Thermodynamics of a charged scalar field in disordered media 
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I. INTRODUCTION AND MOTIVATION 

Disorder plays an important role in the critical behav- 
ior of second order phase transitions [l[ ■ The relevance 
of disorder in the criticality can be assessed qualitatively 
using the critical exponent of the specific heat a for the 
disorder- free system 0; when a > (the specific heat 
diverges at the critical point), the critical behavior of the 
disordered system is changed, when a < (the specific 
heat is finite), disorder has no effect on the critical be- 
havior. On the other hand, at low temperatures quantum 
fluctuations may compete with the random fluctuations. 
For example, the spin-glass ordered ground state - a dis- 
order strongly correlated system - might be destroyed by 
quantum fluctuations. Recently, Carleo et al. [1| demon- 
strated the existence of Bose-Einstein condensation in 
quantum glasses due to frustration. In the present paper 
we focus on the interplay between quantum and random 
fluctuations in a relativistic charged scalar field theory 
with a finite chemical potential. Disorder in relativis- 
tic Bose-Einstein condensation seems not to have been 
considered in the literature, contrary to the case of non- 
rclativistic case, where it has been under intensive study 
since early seminal works 

Disorder has a decisive influence on the zero- 
temperature phase diagram of non-relativistic Bose sys- 
tems. As emphasized by the literature, there is a quan- 
tum phase transition for such systems from an Mott insu- 
lating to a conducting phase. Since no pure Bose system 
can be a normal conducting fluid at zero temperature, 
the conducting-insulator transition must correspond to 
the onset of superfluidity. As shown in Ref. Q, this sce- 
nario is changed dramatically in the presence of a ran- 
dom potential. For the case of a Gaussian colored noise, 
a Bose glass phase also arises and the transition to super- 
fluidity only occurs from this third phase, never directly 
from the Mott insulator. The introduction of a random 
potential in such systems may also imply the destruction 
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of the superfluidity phase, as discussed in Refs. 0,0). In 
particular, a recent study by Lopatin and Vinokur [8|, 
employing the replica method, found a negative shift in 
the condensation temperature of a dilute Bose gas due to 
static disorder - see also Refs. [9l. Hfj|. 

There is an extensive literature on relativistic Bose- 
Einstein condensation (RBEC) following the pioneering 
works of Refs. (TT| - [l6| . which discussed RBEC in flat 
space-times, and Refs. fl~7rl20| . which discussed RBEC 
in curved space-times. While relativistic Bose-Einstein 
condensates are not yet realizable in controllable exper- 
iments like their non-relativistic counterparts, they do 
relate to observable and experimentally accessible phe- 
nomena. One example, of immense current interest, con- 
cerns the condensation dynamics in relativistic quantum 
field theories where creation and annihilation of parti- 
cles play crucial role, like in far-from-equilibrium stages 
of the early Universe and in experiments with relativistic 
heavy- ion collisions [2l|. There is also the possibility of 
Bose-Einstein condensation of pions and kaons |23 - |25| in 
neutron stars. The condensation of these mesons will af- 
fect the equation of state of matter in the interior of the 
star, which has direct consequences on the observable 
mass-radius relation of the star, and will also impact the 
early evolution of the neutron star. 

In real physical situations, the presence of some sort 
of disorder in the system is unavoidable. The disorder 
can be due to uncontrollable disturbances external to the 
system, or due to an incomplete treatment of degrees of 
freedom associated with fields that couple to the field 
of interest. As with nonrelativistic Bose-Einstein con- 
densates of condensed matter physics, one expects that 
disorder will impact the critical behavior of relativistic 
Bose-Einstein condensation. The present study is a first 
step toward a systematic study of disorder in relativistic 
quantum field theory models, in that we focus on a non- 
interacting charged scalar field at finite temperature in 
the presence of nonstatic randomness (the condition for 
staticity will be properly defined shortly). Our model is a 
kind of generalization of the scalar Landau-Ginzburg the- 
ory, where the quenched disorder is described by random 
fluctuations of the effective transition temperature [l|. 

The organization of this paper is as follows. In Sec- 
tion II we present our model. The disorder field couples 
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to charged scalar field via the mass term of scalar field, 
just as in the random-temperature Landau-Ginzburg 
model. We consider weak disorder and implement a per- 
turbative expansion for the free energy as power series 
expansion in the strength of the disorder field. In Sec- 
tion III we study the thermodynamics properties of the 
free relativistic Bose gas at finite density with random- 
ness. In Section IV we obtain the critical temperature 
in the presence of random fluctuations. Conclusions and 
Perspectives are presented in Section V. The paper in- 
cludes an Appendix containing details of renormalization 
of the propagators. Throughout the paper we employ 
units with h= c = ks = 1- 



II. SCALAR FIELD AT FINITE 
TEMPERATURE IN DISORDERED MEDIA 

We are interested in studying the effects of randomness 
on a charged scalar field ip of mass m in equilibrium with 
a thermal reservoir at temperature T. We employ the 
imaginary time formalism of Matsubara [2(| to write the 
partition function of the model in the grand canonical 
ensemble as [U, H3 



Z=[N(/3)] 2 / [D<p][D(p*] e s[v - v '\ 



(1) 



where the action S[ip, ip*] reads 



S[p, ip*] = dr dx 
Jo Jv 



(d t + ifi) ip* (d t - ip) ip 



-\7<p*X7<p — m 2 ip*ip — \(ip*ipY 



(2) 



where V is the volume of the system, /3 = 1/T, p, the 
chemical potential associated with the conserved charge, 
and dt = id T . The field ip satisfies the Kubo-Martin- 
Schwinger [28| boundary condition tp(r, x) = c/j(t + /5,x). 
iV(/3) is a /3-dependent but /^-independent constant that 
comes from the integration over the canonical momentum 
conjugated to the field ip [27|. In this first study, we 
consider the non-interacting case A = 0. 

Next we consider the effects of a random noise 
source over the quantum matter field as in the random- 
temperature Landau-Ginzburg model, but generalized to 
the case of r-dependent noise. That is, we perform the 
replacement m — > m 2 (l + v), where v = v(t, x). The 
partition function given by Eq. ([!]) becomes replaced by 



Z[v] = [N(f3)f / [D<p] [Dip*] e^Thv.v*^ 



where 



S T [v, ifi, ip*] = S[ip, ip*] + Si[u, ip, ip*], 



(3) 



(4) 



with S[ip,ip*] given by Eq. @ with A = and the con- 
tribution due to disorder is given by a random-potential 



term 

Sj[v, p, ip*] — — m 2 I dr I dx. u(T,x)ip*(T,x)ip(r,x). 
Jo Jv 

(5) 

The physical picture is that the random fluctuations de- 
scribe average effects of external disturbances on the sys- 
tem or of degrees of freedom associated with fields that 
have been integrated out formally. Although similar to 
time-dependent noise (29j . a r dependence as in above 
should be understood as arising naturally when inte- 
grating out fields in favor of effective interactions of ip. 
Thence, we are not considering a perfect correlation of 
the disorder in the r direction. This is in contrast to 
the anisotropic case where randomness depends only on 
spatial variables. In addition, hereafter by "static noise" 
we mean that one is considering a set of noises such that 
random fluctuations do not depend on r. In this way, 
in the case studied in this paper one may qualify the 
r-dependent noise as being nonstatic. 

The standard procedure to study Bosc-Einstcin con- 
densation is to separate from ip the constant zero mode 



(6) 



where % is a complex field with no zero mode. The x 
field is written in terms of real and imaginary parts as 



1 



(7) 



so that the action can be written as 

Sr[f,Xi,X2,£] =-pVU(0 + S [xi,X2] 

+ Si[v, X i,X2,t], (8) 
with the classical potential 

U(0 = (m 2 -V 2 )e, (9) 
and the quadratic part 

So[xi,X2] = ~ J dr J dx ^d T xid T xi + VxiVxi 

i / 2 2\ 2 i / 2 2\ 2 

+ (m - n )xi + (m - /i )x 2 

+ d T X2drX2 + Vx 2 V%2 



2«m(X2<9 t Xi - Xi<9tX2) 



(10) 



In Eq. (|10p we neglected the linear terms in the field x, 
because their contributions will be proportional to terms 
like x(p = 0) = 0. In turn, the random contribution is 
given by 

Si[v,Xi,X2,i]=-m 2 J dr J dx|^ 2 ^(r,x) 
1 



^ v(t, x) |^xi (t, x) + x 2 (r, x) 
v / 2^(t,x)xi(t,x)J. 



(11) 
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Here we consider the random function v(t, x) as a 
Gaussian distribution given by 



P["]=Po cx p(-^/ 



d d x{v{x)f 



(12) 



where x = (r, x) and po is the normalization constant of 
the distribution. The quantity a 2 is a parameter associ- 
ated with the intensity of the disorder. We wil l denote 
the mean value over the random variable as (• • • ), defined 
by 



A[ V ]= [Du]P[v]A[v], 



(13) 



with A[v\ being any functional of v. From Eq. (|T2"j) . we 
have a white noise with two-point correlation function 
given by 



u(t, x>(t' , x') = a 2 S(t - t')S 3 (x - x') 



(14) 



As well known, it follows from the Gaussian distribution 
that 



v(xi) ■ ■ ■ v(x 2n +i) = 0, 



(15) 



v{xi) ■ ■ ■ v(x 2n ) = £ \ v{x )i>(x k ), (16) 

pair comb, pairs 

where n is an integer. 

We are interested in calculating the critical tempera- 
ture of condensation. We follow closely the path used for 
the noiseless case [3(3|, in that the critical temperature 
is determined for a fixed charge density while simultane- 
ously fixing the variational parameter £. Noise average 
is taken into account using Eq. (fT3|) . with A being the 
charge density p = Q/V. Specifically, the density is the 
derivative of the free-energy with respect to p, and there- 
fore: 



[Dv\P\v\p = [Du]P[u] 



1 d\nZ R [v] 
W dp 



_t__d_ 

Wdp 



f[Du]P[v] ]nZ R [u] 



1 d]nZsM (m 
/3V dp ' { ' 

Note that we are considering a situation where one has to 
deal with two kinds of averages, namely thermal averages 
and noise averages, which are not treated on the same 
footing. This can be justified when the characteristic 
time scale of the change in disorder is much larger then 
the time of observation of phenomena of interest. 

Equation (fT7|) requires a method to evaluate the aver- 
age over noise realizations of the free energy. For static 
noise and arbitrary noise intensities the replica-trick is 
widely used [l|. Here we consider the weak- noise limit 
and use a perturbative approach [3ll |32| , in that one ex- 
pands the partition function in a power series in the noise 
v. This will be discussed in the next Section. 



III. THE PARTITION FUNCTION AND NOISE 
AVERAGE OF THE FREE ENERGY 

In the weak-disorder limit, the partition function in 
Eq. ([3]) can be expanded in a power series in Sj\ 

Z[v] = (N{P)f / [d X i] {d X2 }e-P vu ^ s ° £ ^, (18) 



where 5/ = Si[u, xi, X2, £] and S Q = 5 [xi,X2]- Taking 
the logarithm of both sides we get 

In Z[v\ = -PV U(0 + In Z + In Z t \v\, (19) 

where the disorder-free part is given by 

(N(f3)) 2 [[D Xl ][D X 2]e s A, (20) 



(21) 



In Zo = In 

and the corrections due to disorder are given by 



l^M=ln(l + £^). 



Here the averages ((• ■ ■ )) arc defined using the unper- 
turbed ensemble represented by the action Sq: 



<(•••)> = 



j[D Xl ][D X 2]e s ° 



(22) 



We now proceed to evaluate Eq. (JT9J) up to second 
order. We start with the zeroth-order term, given by 
Eq. (I2TJ1) . This quantity is calculated explicitly in Ap- 
pendix [A] and the result is 



lnZ = - V 



dp 



(2tt) 3 

+ ln (l_ e -/3Mp)+/^ 



/3w + ln(l-e-' 3Mp) - Al) ) 



(23) 



where w(p) = \/p 2 + m 2 . Next, expanding Eq. (1211) up 
to second order in the noise field, one obtains 

\nZM = {Si) + \[{Sj) - {S,) 2 ^ (24) 

where Si is given by Eq. (fTTjl . Let us now consider the 
stochastic average of such a contribution. From Eq. (|15p , 



we have that (Si) = 0. The other terms are obtained 
using Eqs. dTTJ) and (fT4"|) : 

LTZTR = mV 2 j" dr ^ dx e { X \) + 1 (< X f) 

+ (xt) - (xi) 2 - (xi) 2 



+ 2( X l X D-2(x\)(xl)) 



(25) 
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The derivation of the ensemble averages in Eq. ([25]) is 
outlined in Appendix [Bj and the result is 



\slZi\v] = (0V)m 4 a 2 
+ 1 



dp 

(2^ 



25?l(Wn,p) 



+ 



1 fly 



dp 
dp 

(27T) 5 



2>?iK,p) 



£>22 (^n,P) 



(26) 



where 2?^(w„,p), i,j = 1,2 are the zero-order propaga- 
tors of the fields x j ■ Since the propagators have divergent 
vacuum contributions, Eq. (|26p must be carefully regular- 
ized. The renormalization of the propagators is discussed 
in Appendix [C] After carrying out such a procedure we 
get 



lnZ/M = (pvy 



n m (2£ 2 + n m )-n 2 , (27) 



where the quantities Tl v and TL m (0, p) arc obtained in 
Appendix [Cl they are given by 

dp 



n„ = - 
2 



and 

n m = n m (j0,/i) 
1 

+ 



1 

(27T) 3 w( P ) 
dp 1 

(2^) 3 w(p) 



(28) 



1 



c £Mp)-m) 



,/3(w(p)+/x) _ x 



The last term in Eq. ([27)1 is a divergent contribution to 
the zero-point energy density. Since in non-gravitational 
physics only energy differences are measurable, one can 
normalize the vacuum energy density at zero and then 
the renormalizcd In Z up to second order in the random 
noise will be given by: 



1 



— lnZ R [u] = [p 2 - 



dp 



(2tt)= 



In 1 



i 2 +m 4 cr 2 U m (B,p)] e 



In 1 - 



-/3(w(p)+ M ) 



U 2 m (3,p). (30) 



Notice that thermal contributions arising from the cor- 
rections due to random fluctuations have a similar quali- 
tative behavior for large (or small) as the unperturbed 
part. In addition, notice that the momenta integrals are 
convergent for \p\ < m. In the next Section we discuss 
the determination of the critical temperature. 



IV. THE CRITICAL TEMPERATURE 

As discussed in Ref. [3(|, the parameter £ is not de- 
termined a priori and it should be treated as a varia- 
tional parameter, related to the charge carried by the 



condensed particles. At fixed 3 and p, the free energy is 
an extremum with respect to variations of £. The deriva- 
tive of the above equation with respect to £ implies that 
£; = unless 

p 2 - in 2 + m 4 a 2 II TO (/3, p) = 0. 

Considering a power series expansion for p in the effective 
coupling m 4 a 2 we have, at first order 



\H\ = [m 2 -mVn ro (/3,m)] 



1/2 



or, since we are working in the weak-disorder limit: 



1 



2nV 



II m (/3,m) 



(31) 



(32) 



However, like in the noiseless case, this variational condi- 
tion does not determine the parameter £. In order to fix 
its value the appropriate route is to calculate the charge 
density p = Q/V, which will also enable the determina- 
tion of the critical temperature. As discussed previously, 
the charge density is given by 



P = 



1 

W 



dlnZ R [ 
dp 



fi—m' 



where 



1 - 



2m' 



IL m (0,m) 



(33) 



(34) 



(29) Using Eq. ([50|) . the above expression becomes: 



P = 



2m 



m 4 a 2 



A(Am') 



e + P*W,m') + Pl (l3,m'), 
(35) 



where the unperturbed thermal contribution is 



p*{p,m') = 



dp 

(2tt) 3 
1 



1 



e /3(w(p)-m') _ i 



and the quantity pi{3,m') is defined as 
m 4 o~ 2 3 

Pi(P,m') = —-^A(8,m')IL m (8,m'), 



(36) 



(37) 



with 

A(/W) 



dp 1 

^(p) 
1 



1 



sinh 2 (^(w(p) - m')/2) 
(38) 



sinli (/?(w(p) +m')/2) 

If the density ~p is held fixed and the temperature is low- 
ered, p will increase until the point p = m! is reached. 
If the temperature is lowered even further, then p* + pj 
will be less than p and so 



^ =4) p-[p*(0,m')+p I (0,m')] 



8m' + m 4 cr 2 /3A(/3,m') 



(39) 
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Hence, the critical temperature j3 c 1 is determined im- 
plicitly by the equation given by 



p = p*(Pc, m') + piiPc, m). 



(40) 



In order to clarify the influence of random fluctua- 
tions over thermal fluctuations, let us study the behavior 
of the critical temperature under two different circum- 
stances. First, let us analyze the nonrelativistic limit, 
(/3 m 1) of Eq. (f4"0"l) . In this limit the contributions 
from the antiparticlcs are negligible and, since p <C m, 
we get w(p) — m+p 2 /2m and w _1 (p) ~ m _1 . Therefore 
we get 



n m (/3,m') = -^ ( 
2m J 



1 



dp 

J2nf e /3[£(p)+7] - 1 ' 



and 



dp 



(41) 



(42) 



1 



sinh 2 {/3[£(p)+7]/2} 



4 d 



II OT (/3,m,7), 



where E(p) = p 2 /2m and m' = m — 7, with 



/3 7 = /3- 



2m 



II m (/3,m) 



m 4 g 2 C(3/2) 
4(/?m)V2 (2tt) 3 / 2 



(43) 



(44) 



We performing the above integrals employing spherical 
coordinates, and using result [33J 



1 xp- 1 r(p) 
dx ~^ — = — ^ L w)> 



(45) 



where T(x) is the usual Gamma function and Li„(z) is the 
polylogarithm function, whose series definition is given 
by 



(46) 



Therefore, inserting the results in Eq. (|4U|) one has for 
the total charge density 



2tt/3 c 



3/2 



Li 3/2 (e-^ 7 ) 



8tt 2 



2/3 c m 



1/2 



Lii/2 (e 



(47) 



Since we are working in the weak-disorder limit, we may 
use the following power series expansion about z = 0, 
valid for all s ^ 1,2,3,- •• [H: 



Li s (e- z )-r(l-s)z s - 1 + ^ 



M) fc C(s-fc) „ fc 



fc! 



z fe , (48) 



where £(x) is the usual Ricmann zeta function. There- 
fore, keeping only terms up to order m a 2 : 



2tt/3 c 



3/2 



C(3/2) 



1 

1 - - 77 
2 ' 



1/2 



1 



where 



77 = r){/3 c ,m) 



(/3 c m)V2 (2^)3/ 2 C(3/2)' 
Hence, the critical temperature reads 



p 



2?r 



m VC(3/2) 



2/3 



1 1 1/2 17 
1 + 3^ +36^ 



(49) 



(50) 



(51) 



Equation (|5Tj) is the expression for the critical tempera- 
ture in the nonrelativistic limit. Solving such an equation 
by iteration one arrives at a power series expansion of /3~ 1 
in the effective coupling m 4 er 2 . Therefore, one has 

Pc 1 = Pn 1 {l + I Wn.m)] 1 ' 2 + g 77(/3 n ,m) J , (52) 

where the nonrelativistic critical temperature f3~ l in the 
absence of disorder reads 



2tt 



m VC(3/2) 



2/3 



(53) 



We see that the effect of disorder in the non-relativistic 
limit implies in a higher condensation temperature than 
the one in absence of disorder. 

Now let us analyze the ultrarelativistic limit (Jim <C 1) 
of Eq. ([4*0")) . In this limit we have p ^> m, so that one has 
Lj(p) ~ p. Then 



p*((3,m') w h(m') - h(-m'), 



(54) 



and 



Pi(P,m') 



m A a 2 (3 



F(m') + F(-m') 



16 

G(m') - G(-m') 



where 



M m ') = 7^3 / 



e (3(p-rn') 



2 Z" 00 r,pP(p-™-') 
G(m') = 4 / dp 



(55) 



(56) 



(57) 



77 Jo [ e /3(p-™') - 1] 

/3 dm 1 ^ ' 



(58) 
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and 



1 - 
m < 1 — 



Am 2 



II m (/3, m) 
F(m) + F(-m) 



(59) 



First let us consider F(m'). We use Eq. (|45|) and the 
result (H: 



Li n (e 2 ™ z ) + (-!)" U n (e' Amz ) 



{2m) n 



B n {z), (60) 



where B n (z) are the Bernoulli polynomials, which have 
the following series representations [35| 



, . (-l)"2(2n- 1)! ^ S in2fc7ra; 
B 2n _i(x) = /0 _, 2 „_ 1 2^ 



2n-l ' 



and 



Hence 



( 27r )2n-i ^ k 

x ' k=l 



„ , , (-1)™- 1 2(2n)! ^ cos2fc7rx 

S 2n(x) = 2^ ■ 



(27r) 2n 



fc=i 



F(m') + F(-m') = -^B 2 (x), 



(61) 



(62) 



(63) 



where 2itix = /3m'. In the same way we obtain for G(m' 

Ai 



G(m')-G(-m') = -^»B 1 (x). 



(64) 



Finally, following similar steps as above we obtain for 
h{m'): 



4-7T7 

h(m')-h(-m') = ^B 3 (x) 



(65) 



Therefore, inserting the results in Eq. (|40[) . one obtains 
for the total charge density: 



P 



Airi 
3/3^ 



■Bn 



/3m' 
~2~xi 



Aizfr 



2m 



2m 



(66) 

A more transparent expression can be obtained when we 
make use of the relations in Eqs. (j6"T|) and (f6"2"j) for the 
Bernoulli polynomials. Since /3m <C 1 and we are work- 
ing in the weak-disorder limit, one may Taylor expand 
the circular functions in such equations and neglect all 
terms of order (/3m') 2 and higher. Inserting the results 
in Eq. ([66)) we get, up to order m 4 <7 2 : 



P 



i 



16tt 2 



24(/3 c m) 2 



(67) 



Hence the ultrarelativistic critical temperature in the 
weak-disorder limit is given by: 



^ = 



32tt 2 48(/3 c m) 2 



(68) 



This is the result for the critical temperature in the ul- 
trarelativistic limit. As above, solving such an equation 
by iteration one arrives at a power series expansion of 
fi^ 1 in the effective coupling m 4 cr 2 . Therefore, at first 
order, one has 



Pc 1 = Pu 1 



32tt 2 



48(/3 u m) 2 



(69) 



where the ultrarelativistic critical temperature j3 u 1 in the 
absence of disorder is given by 



(B 



1/2 



(70) 



As in the nonrclativistic limit, we find a positive shift in 
the critical temperature due to random fluctuations. 



V. CONCLUSIONS AND PERSPECTIVES 

In this paper we investigated the effect of weak disorder 
on a free relativistic charged scalar field in thermal equi- 
librium with a reservoir. We studied the effect of coupling 
of a random field to the scalar field in the situation where 
Bose-Einstein condensation takes place. We considered 
a quenched disorder which couples linearly to the mass 
of the scalar field, just as in the random-temperature 
Landau-Ginzburg model. After performing noise aver- 
ages of the free energy, we obtained the corrections to the 
critical temperature for the free Bose gas at finite density. 
We showed that the effect of the randomness is to increase 
the critical temperature in both nonrelativistic and ul- 
trarelativistic limits. Naturally, one should keep in mind 
that this is a result valid for weak disorder and obtained 
in the framework of a perturbative expansion in the noise 
intensity. It remains to be seen if the same result is at- 
tainable with a nonperturbative calculation, e.g. using 



a replica trick. For static (r-independent) noise, appli- 
cation of the replica-trick consists in the following 
Using the fact that InZ = lim„,_>o (Z n — l)/n, one has 
that ]nZ[u] = \im n ^o (Z n — l)/n, where Z n = Z n \y\. 
The Z n is interpreted as the partition function of a new 
system, formed from n statistically independent copies of 
the original system. The quenched free energy functional 
is defined as F q (h) = F q = lim n ^o Z n — l/n, showing 
that the quenched free energy functional can be calcu- 
lated from a zero-component field theory. 

To conclude, we remark on an important point with 
respect to the fact that random mass models generate 
effective interactions that mimic a negative coupling con- 
stant. Many authors claim that non-rclativistic bosons 
only make sense in a random potential when they present 
repulsive interactions @. Nevertheless, there are many 
examples that justify the study of a free theory in the 
presence of random fluctuations. For instance, relativis- 
tic scalar field models with negative coupling constant 
were investigated in the literature and meaningful results 
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were obtained - see for example Refs. |36H41|. Based 
on the results obtained in Rcf. [52], where it has been 
shown that the theory with a negative coupling con- 



stant develops a condensate, Arias et al 
the thermodynamics of a asymptotically free Euclidean 
self-interacting scalar field defined in a compact spatial 
region without boundaries. In addition, based on the re- 

1 I I L 

suits of Refs. [4J, |45|, we have discussed the effects of 
randomness over quantum fields in different physical sce- 
narios. Ref. [3l| proposed in condensed matter physics 
setting an analog model for fluctuations of the light cone. 
Also, a free massive scalar field in inhomogeneous ran- 
dom media was studied in [32|]. After performing the 
averages over the random functions, the two- and four- 
point causal Green's function of the model were presented 
up to one-loop approximation. More recently, Refs. (46j 
and [47| investigated the influence of fluctuations of the 
event horizon on the transition rate of a two-level sys- 
tem which interact with a quantum field. The natural 
extension of such works is to consider the renormalized 
vacuum energy density associated with a quantum field 
in the presence of a stochastic ensemble of fluctuating 
geometries. This is under investigation by the authors. 



(Brazilian agencies) 



43 1 discussed Appendix A: Calculation of the partition function 



In this Appendix we calculate the logarithm of the free 
partition function given by Eq. (f^0|) . We start by intro- 
ducing Fourier series to the fields xi and \2- 



X*(t,x) 



v) EE' 



V 



1/2 



,i(p-x+oj„T 



) X l ;n(p), (Al) 



with i = l,2 and j3u> n = 2irn due to the constraint of 
periodicity Xi(0,x) = Xj(/3,x) for all x. Inserting this 
last result in the free field action given by Eq. ([TOf we 
obtain, after performing an integration by parts: 

«,P \ / 

(A2) 

where we discarded a total derivative term and we defined 
the matrix 9 = 0„(p) as 
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(A3) 

with w(p) = \/p 2 + m 2 . Thus, using Eq. ()A2|) the loga- 
rithm of the free partition function now becomes 

lnZ = ln(Af(/3)) 2 +lnJ(/3, A i) (A4) 

where 



J(P,f*) =nil / d X2-,n(p) cx p j~2 22 X 2 ;"( p )x 2 ;-"(~ p ) 



X / dXl;n(p) exp-{ --\Qn(p)\ n Xl;n{p)Xl;-n(-p) 



x exp 



Qn(p) X2;n(p)Xl;-n(-p) - r ©"CP) Xl;n(p)X2;-n (~P) 

12 2 L J 21 



(A5) 



I 

Noting that Xi;-n(—p) = X*(p)j i = 1, 2, as required by where 
the reality of the fields Xi ( T : x ) i the above integrals are 
just generic Gaussian integrals. Therefore 



J{fi,n) = [det0„(p)] 



-1/2 



(A6) 



where we have discarded an overall constant multiplica- 
tive factor. Inserting this last expression in Eq. (|A4[) . we 
get: 

lnZ = HN{p)f - ilndcte„(p), (A7) 



lndet 9„(p) = In I J] \{ f3 2 [(w„ + z,i) 2 + u 2 (p)] 

^ n p 

+in{nn^ 2 [( w »-^) 2+w2 (p)]i' ( Ag ) 



n p 



and, according to Ref. [2j 



(A9) 
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Hence 



dp 



hiZ = ~VY, J ^ln{/3 2 [( W „4 



dp 



ln{/3 2 [( Wll -z/i) 2 + W 2 ]} 



-2Uln/3^ J 



dp 

(2^)3' 



(A10) 



where lu = w(p) and we consider V to be large compared 
to all other physical lengths so we can replace the sum 
over p with an integral. The frequency sums can be 
performed using standard procedures [48| and the result 
is given by Eq. ([23]). 



Appendix B: Calculation of the corrections to the 
free energy 

In this Appendix we evaluate the perturbative correc- 
tions to hiZ given by Eq. (l25j) . For this purpose, let 
us present some basic derivations. We define the finite- 
temperature propagators in position space as pol ] 

P. y (T,x;r',x') = (xi(T,x) X j(T J ,x')),i,j = 1,2{B1) 

Due to translation invariance, the above propagators 
should depend only on x — x' and r — r'. With the help 
of Eq. (|Alj) . their Fourier transforms are 



2> <3 -(w n ,p)= f dr [ dxe-^^^^f^x) 
Jo Jv 

= /3 2 (Xi;n(p)Xr,-n(-p)), (B2) 

where T>ij(r, x) = 2?^ (t, x; 0, 0). Alternatively, we note 
that the propagators can be expressed as functional 
derivatives of the partition function: 



V ij {u n ,p) = 2(V%) U. 



(B3) 



Hence, with the help of Eq. (|A4p the zero-order propa- 
gator Pj^Wnjp) is given by 



5n(u ni p) 



-2f3 2 



S In Zo 

5[e„( P )] 



ii 



(B4) 



Similar expressions hold for the other zero-order propa- 
gators. Using Eq. (|A7|) together with Eq. (|A3p . we get 



1 r 



Vl 1 (w n ,p) = - V + {uj n ,p) + X>_(u;„,p) 



where 



£>±(u; n ,p) 



K±v) 2 + w 2 (p)' 



(B5) 



(B6) 



In the same way: 



and 



£>22 K » P) = 2 P + ^ n ' P) + X> - ( W « , P) 



(B7) 



V\ 2 (w n ,p) = -D^(w„,p) 



K + (w(p) + ^) 2 K + Mp) - m) 5 



(B8) 



We may apply the above technique in order to calculate 
the various ensemble averages emerging in Eq. (|25|) . Us- 
ing Eqs. (|A1[) and (|B2|) one can easily show that 



dr / dx(x?(r,x)>=y;X' 11 (a;„,p). (B9) 

n,p 

Following the same strategy, we have that 

r i 2 



/ dx(xi(r,x)) = — 



_ n , p 



(BIO) 



and 



Y,V° n (u n ,p) 



(Bll) 

Similar expressions are valid for the field %a- Then, using 
the same procedures as above, one has 



£ dr J dx ( X f(r, x))(xl(r, x))= jf " dr jf dx < X 2 (r, x) X 2 (r, x)) = -L 



-1 2 



E^l^'P) 



(B12) 



since Vl x (uj n ,p) = X>S? 2 (w„,p) and En^iaC^p) = 
En ^2i( w n,p) = 0. Inserting the above results in 
Eq. (|25|) and considering the large- volume limit leads us 
to expression (|2"r?|) . This is the contribution to \nZ due 



to randomness up to second order in the noise. Nev- 
ertheless, this is not the final answer; as well known, 
the vacuum contribution to the propagators is divergent. 
Therefore, in order to present a finite result for In Z one 
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must first renormalize ^(wnjp) and 2?2 2 (u; n ,p). This 
we discuss in the next Appendix. 



Appendix C: Renormalization of propagators 

Here we examine the renormalization of the finite- 
temperature propagators 2?$ ) 1 (u; n ,p) and D^KiP)' 
Again following Ref. [3(J we define the self-energy 111 = 
ni(w„,p) with respect to the averaged propagator 2?n 
as 



X>n(w n ,p) = (1+2??! Ei) l V° n 



(CI) 



A similar expression holds for II2 = n.2(w n ,p) which is 
the self-energy with respect to the averaged propagator 
T>22- Hence, recalling Eqs. (|B3j) and (|B4[) . we get, up to 
second order in v. 



(i + ^nO-^i + 22??,^!, (02) 



and 



^l + ^ 2 n 2 )- 1 =l + 2P« 2 ^M. (C3) 



5V. 



22 



Therefore, inserting Eq. (|26[) in the above equations and 
expanding their left-hand sides to first order yields the 
following expressions for the self-energies 



n 1 (u,„,p)=-2mVe 2 



4 2 

ma 



E 



dp 



2?° 1 K,p), (CM) 



and 



n 2 K, P ) = / (^1 ^ 2 K,P). (C5) 



On the other hand 



1 \r- f dp 1 f dp 1 



1 



1 



e /3(u( P )-/j) _ ^ e ^(w(p)+/i) _ 1 



(C6) 



for i = 1, 2, where we have performed the sums over the 
Matsubara frequencies in the usual way [H] . In this way 
we have 



ni(w„,p) = -mV [2£ 2 + it, + n m (/3,/^)] 

and 

n 2 K, p) = -mV 2 [n„ + n m (/?, /x)] . 

where we have defined 



n, 



1 f dp 1 



(27r)3 w (p) ; 



(C7) 
(C8) 

(C9) 



and 



,„ . 1 f dp 1 
n «^)=2j/ (2^^p) 



e ^(^(p)-M) - 1 



3 ^(w( P )+/j) _ ^ 



(CIO) 



Since IT, is a divergent quantity, in order to avoid physi- 
cally meaningless results the following counterterm must 
be added to the original action: 

5S= —8m 2 I dr dx<p(p* 
Jo Jv 

= - Sj f [dr j^^e + xl + xi), (Oil) 



I 

where we have droped terms linear in xi and Xi- Treating 
this as an additional interaction, we see from Eq. (|2"Tj) 
that to lowest order this counterterm contributes to In Zj 
as 



W)5m 2 e J -^ [dr J v d X (( X l) + { X l)). (C12) 



The counterterm should be chosen so that 



Sm 2 - m 4 a 2 n v = 0. 



(C13) 



In this way we get a finite result for the propagators. 
Whence, collecting the above results, the contribution to 
In Z up to second order in the noise will be 



4 2 
ma r 



biZi[v] = (j3V)—j- n m (2^ + n m )-n^ .(ci4) 
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